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1. Introduction 

The vacuum of the quantized Yang-Mills gauge field theory contains most information about 
its distinct features like confinement, chiral symmetry breaking, etc. In the Hamiltonian formula- 
tion, this information is carried by the ground-state wavefunctional. In temporal gauge, in (d + 1) 
dimensions, the problem looks very simple: One strives to solve the Schrodinger equation 

J^^>o[A] = J d''x!^-^^^^^ + \Ff^{xf^^'o[A]=E,^'o[A] (1.1) 

with an additional constraint that physical states are invariant under infinitesimal local gauge trans- 
formations (GauB' law): 

S-'^dk + ge'^'^Ai) -^nM = 0- (1-2) 

Subtleties of this problem were outlined more than 30 years ago by one of us (J.G.) [jl]]. It 
was argued that at large distance scales one expects the wavefunctional to assume the so called 
dimensional-reduction form: 



'¥f[A] « exp 



-fLjd'xFtj{x)Ftj{x) 



(1.3) 



i.e. a vacuum with color-magnetic fields fluctuating independently in each spacetime point. In such 
a case the computation of a spacelike loop in (t/ + 1 ) dimensions reduces to the calculation of a 
Wilson loop in Yang-Mills theory in d Euclidean dimensions. If the property existed for Yang- 
Mills theories (3 + 1 ) and (2 + 1 ) dimensions, then these would be confining, since the theory in 2 
euclidean dimensions exhibits the area law. ^ However, the dimensional-reduction form cannot be 
the whole story, it does not provide correct short-distance structure of the theory. 

As a step forward, Greensite proposed a systematic strong-coupling expansion of the Yang- 
Mills vacuum wavefunctional in the form: 

»Fo[t/]=^exp(/?[[/]), (1.4) 

where the function R in the exponential is an expansion in terms of closed loops - products of link 
variables U along closed contours on the lattice: 

R[U] = 



+ larger contours. 

(1.5) 

It was later shown by Guo, Chen and Li ^ that for smoothly varying gauge fields the first terms 
of the expansion are expressed through the color magnetic field strength, B"{x) = F"2{x), and the 
covariant laplacian in the adjoint representation, S''^ = &k- ^k, where S>k[A] denotes the covariant 
derivative in the adjoint representation: 

R[U\ oc -1 (aK-oTr [B^]-a^K2Tr [B{-&^)B] + . . .) , (1.6) 



'in the rest of this paper, we will discuss exclusively the case of (2+1) dimensions. 
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where 

Kb = ico + 2(ci+C2 + C3), K2 = {ci with Cq = ^(jS^), d , C2, C3 = ^(jS^) . (1.7) 

With only a few exceptions (see references in [Q]), there was not much work done in this 
area after the initial efforts. However, recently interest in the problem of the Yang-Mills vacuum 
wavefunctional has been revived and various plausible proposals for the vacuum states have been 
advanced. We will present tests of two of them which differ by their motivation but share some 
common features. More results and other proposals will be covered elsewhere [Q]. 



2. Two proposals 



With a grain of imagination one can assume that an expansion of the form ( |1.6| ) might result 
from a vacuum wavefunctional 



»Fo[A]=^exp 



J cfxcfyB\x) K'^[-&'^] B\y) 



(2.1) 



with a kernel depending on the gauge coupling g and the adjoint covariant laplacian. That kernel 
cannot be arbitrary - one should be able to reproduce the QED vacuum wavefunctional [for SU(2) 
in (2 + 1) -dimensional spacetime] for g — 0, which is known to be 



^o[A] oc exp I - ^ j (fxcfy [V x A(x)] {-jl=^ [V x A 



iy)] 



The kernel has thus to satisfy the condition: 



\imK[ 



ah r 



(2.2) 



(2.3) 



The above condition is, by construction, satisfied by the approximate vacuum wavefunctional 
proposed recently by the present authors : 



*PGo[A]=^exp 



i j (fxd^yB"{x)\ 



ab 



(2.4) 



Here Ao denotes the lowest eigenvalue of (— ^^), and m is a constant (mass) parameter proportional 
to ~ l/j3. It is similar to a proposal by Samuel [^; the difference lies in the subtraction of Xq 
which is crucial, since (— has a positive-definite spectrum and we have hints that its lowest 
eigenvalue diverges in the continuum limit. 

The supporting evidence for this proposal comes from four sources [^, ||]: 

1. The proposed form is a good approximation to the true vacuum for strong fields constant in 
space and varying only in time. 

2. If we divide the magnetic field strength B{x) into "fast" and "slow" components, the part 
of the vacuum wavefunctional that depends on Bsiow takes on the dimensional-reduction form. The 
fundamental string tension is then easily computed as 



Of = 3mg 1 16. 



(2.5) 



3 



In quest of the Yang-Mills vacuum wavefunctional 



Stefan Olejnik 



3. If one takes the mass m in the wavefunctional as a free variational parameter and computes 
(approximately) the expectation value of the Yang-Mills hamiltonian, one finds that a non-zero 
(finite) value of m is energetically preferred. 

4. Results for the mass gap, and the Coulomb-gauge ghost propagator and the color-Coulomb 
potential computed from *Fgo [A] are in good agreement with those derived from standard Monte 
Carlo simulations (see below). 

Despite this evidence, our proposal represents only an educated guess, motivated by the form 
of the QED vacuum state, dimensional reduction, and gauge invariance. 

A more sophisticated approach has been followed in (2 + 1 ) dimensions by Karabali, Nair and 
collaborators [^. In the temporal gauge (Aq = 0) they combine the remaining two components of 
the gauge potential into complex fields: 



A = i(Ai+/A2), A = i(Ai-/A2), (2.6) 
and then introduce new variables, a matrix-valued field M G Sh{N,'if), which is related to A, A via 



A = -{d,M)M-\ A = Mt \d^^^), 



(2.7) 



where z = xi — ixj and z = xi + ixj are the usual holomorphic variables in the complex plane. 

Under a gauge transformation Q., M transforms covariantly, M — IIM, and can be used to 
define gauge-invariant field variables: 



H = MTM, 7"~Tr(r(5H)H-^) 



(2.8) 



through which one can express the hamiltonian, inner products of physical states, and the vacuum 
wavefunctional. 

Karabali et al. argue that the part bilinear in field variables J", when expressed in usual vari- 
ables, has the form: 

ab 



*Pkkn[A] « exp 



d^xd^y B"{x) 



1 



\/-V2 + + 



m 



B\y) 



xy 



(2.9) 



which, however, is not gauge-invariant. One can imagine that higher-order terms in J" might 
convert the ordinary laplacian in Eq. (2.9) into the covaiiant laplacian, leading to: 

ab 



*I'kkn'[^] ~exp 



d^xdyB'\x) 



1 



m 



B\y) 



(2.10) 



This form is still hardly sustainable, because of the divergence of the lowest eigenvalue Xq of the 
adjoint covariant laplacian (discussed above). 

Instead of the form ( 2.10| ) we will subject to lattice tests a "KKN-inspired" or "hybrid" wave- 
functional which has (2^) and ( 2.1C ) as starting point and agrees with them for abelian gauge 
configurations, but in which, similarly to the GO proposal, the covariant laplacian (— is re- 
placed by the subtracted one (— — Ao): 



^hybrid [A] =^ exp 



\ j d^xd^yB"{x)\ 



1 



ab 



^/{-^'^-'ko)+nfi + m 



xy 



(2.11) 
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3. Tests of the proposals 



Our aim is to test how good/bad is the approximation of the true Yang-Mills vacuum wave- 



functional by the proposed approximate forms, Eq. ( [2.4[ ) and ( |2.1 1| ). To achieve this goal, we take a 
set of operators {2 [A]} that depend on gauge fields A, and compute (and compare) their expectation 
values: 

• (*Fq™''|2[A]|*Fq™'^) = {Q[A] )mc in Monte Carlo lattices, i.e. an ensemble of two-dimensional 
slices of configurations generated by MC simulations of the three-dimensional eucUdean SU(2) 
lattice gauge theory with standard Wilson action at a coupling jSe; from each configuration, only 
one (random) slice at fixed euclidean time is taken; 

• (*I'o|2[A]|*I'o) = {QiM )recursion in "recursiott" lattices, i.e. an ensemble of independent two- 
dimensional lattice configurations generated with the probability distribution given by a proposed 
vacuum wavefunctional, with parameters m and fixed to some reasonable values, to be able to 
compare to the MC ensemble. 



Numerical simulation of |*Fop The generation of recursion lattices whose probability distribu- 



tion P[A] is given by the square of a wavefunctional of the type (2J) is based on the following 
idea [^] : Define a probability distribution for gauge fields A with the kernel K controlled by an 
independent "background" configuration A' 



^[A;K[A']]=^exp 



J cfxcfyB"{x;A)K"^[A']B\y;A) 



(3.1) 



where the field strength B is computed from A, and both A and A' are fixed to an appropriate gauge. 
If the variance of the kernel K[A] in the probability distribution P[A] is small after the choice of 
gauge, then one can write down a chain of approximate relations: 



P[A] = ^[A;K[A]] ^ ^[A; {K)] = ^ 



A; dA' ^[A;K[A']]P[A" 



dA' 3^[A;K[A']]P[A' 



(3.2) 



The probability distribution P[A] can then obtained by solving ( |3^ iteratively: 

p(i)[A] = ^[A;^[0]], P^^+^\a]= j dA' ^[A;K[A']]P^^\a']. (3.3) 

Practical implementation of the recursion procedure consists of the following steps: Choose 
Ai = (axial gauge) and A2 / 0, then 

(i) given A2, set A'2 = A2, 

(ii) ^ [A;/r[A']] is gaussian in B, diagonalize K[A'] and generate a new B-field stochastically, 

(iii) from B calculate A2 in axial gauge and compute everything of interest, 

(iv) go back to step (i), repeat as many times as necessary. 

The procedure converges rapidly, one needs i^(lO) cycles above, and the assumption about a 
small variance of K among configurations is supported a posteriori by the absence of large fluctua- 
tions of the spectrum of K evaluated on individual recursion lattices. 
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Choice of vacuum wavefunctional parameters m and g Our next task is to select appropriate 
values for parameters g (or p = 4/^^) and m of recursion lattices to be able to compare to MC 
lattices with Wilson-action coupling jSe- In our earlier studies ^ with *Fgo we chose j5 = jSe, 
and fixed m using Eq. ( |2.5[ ) by the measured value of the fundamental string tension (in lattice 
units): m(j8E,L) =4jSEaf(i8E,L)/3. 

Another possibility is to use for fixing j8 and m some information about the true Yang-Mills 
vacuum wavefunctional at jSe for a set of simple gauge-field configurations. The square of the 
vacuum wavefunctional for some trial configurations (non-abelian constant fields, abelian or non- 
abelian plane waves) can be computed numerically in simulations of the three-dimensional Yang- 
Mills theory. Take a set of time-independent configurations ^ = {U^^\x),j = l,...,M}. The 
method (proposed long ago [ ]T0| ] and described in more detail in [Q]) is based on the following 
identity: 



\nu'^j^]\^ = y[DU]d{Uo)lls[u{x,0)-U^j\x) 



e-'. (3.4) 



In practice, one measures the probability in a modified lattice Monte Carlo simulation: The links 
at f = are constrained to belong to a configuration from the set In a MC update all links, 
except those with f = 0, are updated by the usual heat bath method. On the t = slice, one of the M 
configurations from the set is selected at random, and then accepted/rejected by the Metropolis 
algorithm. Len Nj denote the total number of times that the j-th configuration from the set is 
accepted, and A'tot the total number of updates of the t = plane. Then: 

|»P[C/W]poc lim ^. (3.5) 

Nto, lai-ge TV tot 

For determining vacuum wavefunctional parameters we measured probabilities of abelian 
plane waves with fixed (maximal) wavelength A = L and varying amplitudes: 



Ui^\ni,n2) = jl-aj{n2)^h + iaj{n2)a3, U^^\ni,n2) = h, (3.6) 



with aj{n2) = J —^cos—^, p=2n-cos — j. (3.7) 

The probabilities measured in the Monte Carlo simulation described above can be parametrized by 
\'i'Mc[U^H^ = expi-RMc[U^j^]-Ro), Rmc[U^'^] = 2{a + yj) (OMcip) + const. (3.8) 
Similarly, for a theoretical Ansatz of vacuum wavefunctional: 

/?Ansatz[f/^^'^] = 2{a + jj) (OAns^M + COnSt. (3.9) 

For the proposals discussed in Section ^:^ 

Okio{p) = — , , ftJKKN(p) = — ; • (3-10) 

Fig. |l] shows results for (Omc vs. p^ in physical units for a number of jSe values and lattice 
sizes L. The scale was set by the conventional value of (0.44 GeV)^ for the physical string tension. 



For abelian configurations, the KKN (2.9) and hybrid ( |2. 1 1| ) wavefunctionals coincide. 



6 



In quest of the Yang-Mills vacuum wavefunctional 



Stefan Olejnik 



3 





fc=6 






0.9 


Pe=^ *- -' 






|5e=12 






0.8 


GO fit 








KKN fit 






0.7 








0.6 








0.5 
















0.4 
















0.3 
















0.2 










y 






0.1 


1 
















0.5 1 1.5 2 2.5 

(GeV^) 



Figure 1: Cumulative data for Omc vs. in physical units, on lattices of extensions L — 16,24,32,40,48, 
and eucUdean lattice couplings jSe = 6, 9, 12. The curves (hardly distinguishable from one another) represent 
COgo{p) and COkknIp) for fitted parameters mphys and ^pj^^j, given in Table [l|. 



i.e. the lattice spacing is a{pE,L) = \/o'f(j8E,.i-)/(0.44 GeV). The data were fitted by functional 
fornis in Eq. ( 3.10 ), the resulting parameters are summarized in Table |l]. 



Parameters of the proposed wavefunctionals in lattice units were 
fixed using: 

/(j8E,L) =gpi,ys a(/3E,L), m(j3E,L) =mphysa(j8E,L). (3.11) 



variant 


2 

'Wphys gphys 


GO 
KKN 


0.771 1.465 
0.420 1.237 



We will present below results for /3e = 9, the actual parameter values Table 1: Vacuum wavefunc- 
used at this coupling are hsted in Table tional parameters (phys. units). 

Results We focus on two important quantities defined in Coulomb gauge, the Coulomb-gauge 
ghost propagator: 



G{R)={{^[A]-X 
and the color-Coulomb potential: 



\x-y\=R 



(3.12) 



\x-v\=R 



\x-y\=R 



y-&[A] 



(-V^) 



1 

V • &[A] 



\x-y\=R 

(3.13) 

It was argued by Gribov [jT3|] and Zwanziger [14|, that the low-lying spectrum of the Faddeev- 
Popov operator, ^[A] = —V • ^^[A], in Coulomb gauge probes properties of nonabelian gauge 



Ik L 




GO 


hybrid 




j3 m 


P m 


9 32 


0.162 0.367 


1.43 0.283 


8.80 0.154 



Table 2: Values of j3 and m for the GO and hybrid wavefunctionals at jSg = 9, derived from the parameters 
in Table |l]. (The value of ^y <yf{j5E,L) comes from Monte Carlo simulations of the standard Wilson action in 
three Euclidean dimensions JI^.) 

■'The results obtained for the GO vacuum wavefunctional with two variants of fixing its parameters, described in the 
text, almost do not differ lUlll. 
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Figure 2: The Coulomb-gauge ghost propagator. Figure 3: | V(0) | in the individual configurations. 



fields that are crucial for the confinement mechanism. The ghost propagator in Coulomb gauge and 
the color-Coulomb potential are directly related to the inverse of the Faddeev-Popov operator, and 
play a role in various confinement scenarios. In particular, the color-Coulomb potential represents 
an upper bound on the physical potential between a static quark and antiquark [|^]. 

An important point to mention is the equality of the vacuum wavefunctionals in temporal 



and Coulomb gauge (see e.g. []160), when evaluated on gauge fields satisfying the Coulomb gauge 
condition V • A = 0, and which lie in the first Gribov region. Our numerical method described above 
generates configurations in the temporal gauge, these are then transformed to Coulomb gauge, and 
Coulomb-gauge observables are evaluated in the transformed configurations. 

Figure ^ displays the equal-time ghost propagator in Coulomb gauge computed in a standard 
Monte Carlo simulation on a 32^ lattice at jSe = 9, together with results obtained from recursion 
lattices with probability distributions given by and ^I'hybnd' generated using /3 and m values 
listed in Table ^ The agreement is quite perfect, for all three ensembles. 

Figure || is crucial for understanding results for the color- Coulomb potential, which are shown 
in Fig. ^. One can evaluate the potential in each individual lattice configuration. Figure ^ dis- 
plays values of \V{0) \ in MC lattices, GO and hybrid recursion lattices. It is clearly seen that most 
configurations have \V{0)\ in the range between about 2 and 10 (about 80%), but there are rare 
instances of configurations with much higher values. These "exceptional" lattices possess a still 
positive, but very small value of the lowest eigenvalue of the Faddeev-Popov operator, and are 
therefore rather difficult do gauge-fix to Coulomb gauge. If we discard from our ensembles con- 
figurations with \V{0)\ greater than some cut K, we obtain results illustrated in Fig. |[ For K" = 5 
or even 10 the agreement of potentials for MC, GO, and hybrid lattices is reasonable. However, 
as the cut is increased, the agreement deteriorates. The GO and hybrid potentials are still roughly 
linear (and hardly distinguishable from each other), but deviate quantitatively from the MC result. 
This indicates a discrepancy in the tails of the probability distributions. While the ghost propagator 



(3.12), containing only a single factor of the inverse Faddeev-Popov operator, is rather insensitive 
to the tails and its values are mainly determined by the bulk of configurations, the color-Coulomb 



potential ( [3.13| ) involves two factors and is more sensitive to the tails. 

But what makes the probability distributions corresponding to *Fqq and ^I'hybnd close? We 
believe that the reason is that both wavefunctionals have - for optimal choices of their parameters 
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Figure 4: The color-Coulomb potential for illustrative values of the cut K. 

- about the same dimensional-reduction limit, and the results are mainly sensitive to that limit, not 
to the detailed functional form of the kernel K that enters the wavefunctional. 



4. Summary 

1. We described a recursion procedure that allows to generate lattice configurations with prob- 
ability distributions given by approximate Yang-Mills vacuum wavefunctionals of certain simple 
forms Uke and ( |2.11| ). 

2. Relative magnitudes of the true vacuum wavefunctional on particular sets of configurations 
(abelian plane waves, non-abelian constant configurations) can be computed numerically. 

3. Parameters of approximate vacuum wavefunctionals can be fixed e.g. by fitting the results 
for long-wavelength abelian plane waves. 

4. The two tested proposals (GO, KKN-inspired hybrid) provide Coulomb-gauge quantities 
almost indistinguishable, and in reasonable agreement with lattice Monte Carlo results (with some 
discrepancy in color-Coulomb potentials). 

5. GO and hybrid vacuum wavefunctionals seem to agree with the true Yang-Mills vacuum 
wavefunctional for the bulk of the probability distribution. 

6. The important common property of both tested approximate vacuum wavefunctionals ap- 
pears to be their almost identical dimensional-reduction form. 

Only a subset of our recent results was covered in the present contribution. The interested 
reader should consult Refs. [Q, ^ |8[ [iT]] for additional details and more data. 
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